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Abstract 

We consider an array of N Josephson junctions connected in parallel and explore the condition for 
chaotic synchronization. It is found that the outer junctions can be synchronized while they remain 
uncorrelated to the inner ones when an external biasing is applied. The stability of the solution 
is found out for the outer junctions in the synchronization manifold. Symmetry considerations 
lead to a situation wherein the inner junctions can synchronize for certain values of parameter. 
In the presence of a phase difference between the applied fields, all the junctions exhibit phase 
synchronization. It is also found that chaotic motion changes to periodic in the presence of phase 
differences. 
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Due to the application of chaotic synchronization in secure communication 
to brain modehng a great deal of investigation has been done in this field. The 
presence of even a small phase diff^erence between the applied fields was found 
to desynchronize a completely synchronized system. Also the phase difference 
was found to have application in taming chaos in dynamical systems. Recently 
it was observed that the end lasers in an array of three laser system was found 
to synchronize while it remained uncorrelated with the middle laser which orig- 
inally connected the two. In this work we study an array of Josephson junctions 
in the presence a phase difference between the driving fields and its effect on 
synchronization and suppression of chaos. 



I. INTRODUCTION 



Chaos in Josephson junction (JJ) has been studied extensively after its presence was 
demonstrated using numerical simulation When we treat JJs within the Stewart- 
McCumber model, the equation describing the behaviour of JJ is identical to the equation 
for a driven damped pendulum which has been studied theoretically for several routes to 
chaos . Thus JJ becomes an ideal physical system to study chaos. The rf- biased JJs 
find practical importance in the construction of devices like parametric amplifiers, volt- 
age standards, pulse generators, SQUID for detection of very weak magnetic fields, etc. 
0,0, ^- For these devices, it is essential to avoid all types of noise, chaos etc. JJs consisting 
of Superconductor-Insulator-Normal metal-Insulator-Superconductor (SINIS) showing non- 
hysteretic I-V characteristics with high damping has been fabricated for programmable dc- 
voltage standards [3] or ac- voltage standards based on synthesis of calculable wave formsQ]. 

Pecora and Carroll in 1990 reported that synchronization of chaotic systems ^ could 
be achieved, since then different types of synchronization such as complete, generalized 
and phase synchronization of chaotic oscillators have been described theoretically and ob- 



served experimentally 
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ll[ |. Synchronized chaotic oscillations have been found in many 



nonlinear systems like lasers, neural network, etc[l^, [l^. Chaotic synchronization also find 
application in communication. It was demonstrated using Rossler oscillators that during the 
transmission of information about a stimulus through an active array, the stimulus created 
the way to be transmitted by making the chaotic elements to phase synchronize [l^. The 
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stability of synchronous state is analyzed by Lyapunov function method 15| and the mas- 



ter stability approach 



16| . Phase difference between the applied fields plays an important 



role in suppressing chaos and the synchronization of chaotic systems. Duffing oscillator was 



studied for the effect of phase difference on chaotic synchronization [17|. Josephson junction 
has been investigated for both periodic and chaotic synchronization. Coupling between self 
generated Josephson oscillations through a microwave transmission line was found to play 



an important role in collective synchronization of JJ array 



18| . In a system of two JJs in 



parallel, the phase difference between the applied fields was found to bring chaotic motion 
to a periodic one for a large range of parameter values |19]. A parallel array of coupled 
short JJs linked together by inductors has been used to fabricate highly sensitive detectors 
[20I . Although there is extensive work on synchronization of coupled JJs, studies on chaotic 
synchronization of JJs is much less. 

In this work we analyze a parallel array of N-coupled JJs with parameters lying in the 
chaotic regime and study synchronization of the system. The paper is organized as follows. 
In section II we discuss the model for an array of JJs linked in parallel with linking resistor 
Rs in between. Section III contains the study of the synchronization in such an array and 
discuss the stability of the synchronous solution. The effect of phase difference between the 
applied fields on synchronization and its role in suppressing chaos is also discussed. Results 
are summarized in section IV. 



JJ, 



JJ, 



JJ„ 



FIG. 1: Schematic representation of an array of J J linked in parallel with a linking resistor Rg. 1 
and 2 are the driving fields. 
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II. THE MODEL 



The equation of a single Josephson junction represented by the resistively and capacitively 
shunted junction (RCSJ) model can be written by solving Kirchoff 's law as 



hC (fcp h d(f) I I I , i\ 



(1) 



where is the phase difference of the wave function across the junction, Zq cos(co't') is the 
driving rf - field and i'^^ is the dc bias. The junction is characterized by a critical current ic, 
capacitance C and normal resistance R. The coupled JJ considered here consists of a pair of 



2l| . A schematic representation 



such junctions wired in parallel with a linking resistor Rg 
of an array of J J wired in parallel with linking resistors is given in fig [H The equation of 
motion for an array of N coupled current driven JJs can be written in the normalized form 

as 



)i + + sin 01 = idc + io cos(fit) — as 



(2a) 
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(2b) 



W + /50Ar + sin 0Ar = idc + io cos(fit) - Us 



■)N-l 



(2c) 



where i varies from 2 to N-1 and the dimensionless damping parameter (3 is defined as 
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R V 2eic 

The normalized time scale is written as t = uJjit' where ujji = {2eici/hCi)^ . The dc bias 
current i'^^ and the rf amplitude are normalized to the critical current i^i. The actual 
frequency u is re-scaled to 1) = u/uji and the coupling factor is defined as as = (Ri/Rs) (3. 

The Josephson junction is chaotic for the parameter values (3 = 0.3,20 = 1.2, = 0.6 
and idc = 0.3. We fix these parameter values for the numerical simulations. The junctions 
are taken to be identical and for a coupling strength of as = 0.37, the outer junctions 
synchronize while the inner junction remain uncorrelated with the two outer ones. It can be 
seen from Fig[2]J^a) that the outer junctions are synchronized whereas Figj2](b) shows that it 
is uncorrelated with the middle junction for an array of three JJs. 



III. STABILITY ANALYSIS 



In order to perform the stability analysis for the synchronized state of N-coupled Joseph- 
son junctions, we first consider three JJs linked in parallel. In the first order form the three 
identical junctions can be written as 

01 = -01 (3a) 

ipi = - sin (pi + idc + io cos(fit) - [^Iji - 1IJ2] 

4>2 = tp2 (3b) 

^2 = -/3V'2 - sin 02 + Ois VPi + ^/'s - 2V'2] 

03 = V's (3c) 

= -/^V^s - sin 03 + idc + io cos (fit + 9) - Us [^/'s - ^2] 

From eqj3a] and [3c] it can be observed that the outer junctions are identical and symmetric 
with interchange of variables in the absence of a phase difference 9 between the applied fields. 
Hence there exists an identical solution for the outer systems given by 0i = 03 = 0(t) and this 
type of behavior where systems show identical behavior is called complete synchronization. 
Due to asymmetry the middle junction may have different dynamics. The stability of the 
synchronous solution of the outer junctions is analyzed by two methods. 

We define the difference variables 0]'3 = ^i^^ and ip^^ = ^^"'^^ and the approximate 
dynamics transverse to the synchronization manifold is obtained by linearizing the corre- 
sponding subsystem consisting of the outer junctions. The equation may be given as 

= ^ils (4) 
^fa = -/5^r3 - cos 0^3 sin 0^3 - asipi^ 

Linearizing eq. H] we get the approximate dynamics transverse to the synchronization man- 
ifold. In terms of the Jacobian matrix we can rewrite the above equation as 

1 

cos 01 —P — as 

where sin 05^3 ^ 0|f 3 and cos 0^3 ~ cos 0i as 0i ~ 03 in the synchronization manifold. The 
eigen values of the matrix are 
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{as + 13) 



4 cos 01 

1 ± W 1 + 



{as + f3f 



(5) 



The stability of the synchronous state is controlled by the eigen values mi^2 [2^. If mi^2 
are complex conjugates with negative real part, the corresponding solution is stable. In the 
above case the average of the term in the radical is found and it is a complex number with 
real part greater than unity. The real part of the largest eigen value is thus found to be 
negative and hence satisfy the criterion for stability of synchronization. 



As a second test, we follow the method given by Landsman et.al [23|] where the conditional 
Lyapunov exponents are calculated with respect to the perturbation out of the synchroniza- 
tion manifold. Eq. [3] reduces to a set of four equations in the synchronized state as the outer 
junctions may be represented by a single set of equations. In terms of the synchronous solu- 
tions (j)(t) and ipit), we can define variables Acpit) = (f>i(t) — (pit) and A'?/'(t) = ipiit) — ip{t). 
Linearizing transverse to the synchronization manifold, we have 



dt 

and 



(6) 



where i=l,3 and J is the Jacobian matrix evaluated at A0(t) and Aipit). Thus we have 

A0i,3 \ 1 \ ( ^h,^ 

A^i,3 J \l -P-as) \ A^i,3 
A0i^3 and A^/'i^s are the perturbations of the outer oscillators from the synchronous solution 

The Wronskian of the linearized system can be related to the trace of the matrix by the 
Abel's formula 



23| 



W{t) 



exp I / (— «s — f3)dt' 







A^ 
A4> Alp 

where we have dropped the subscripts of the linearized variable. The Wronskian gives the 
phase space dynamics of the system. Taking the natural log of the Wronskian we get 

ln[W{t)] = In \A(j)AtP - AV^A0| = - [ {as + (3)dt, (8) 

which is a monotonically decreasing function of time. The sum of the conditional Lyapunov 
exponents is given as. 



M ^ 

Va,- = lim -ln|det(<l>(A0i,3,A^i,3)(t)| (9) 



where $ is the matrix solution of eqns. [6] and [71 The sum of the conditional Lyapunov 
exponents can be now approximated as 

A1 + A2 ^ + (10) 

The sum of the conditional Lyapunov exponents is negative indicating that the phase space 
of the coupled system shrinks to a trajectory representing the synchronous solution. Thus 
the two methods lead to the same conclusion. 
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FIG. 2: (a) The outer juntions are synchronized (b) Outer junction and middle junction is uncor- 
related. The parameter values are f3 = 0.3, io = 1.2, u; = 0.6, i^c = 0.3, = 0.37. 

Now we analyze the subsystem constituted by the outer and the middle junctions. We 
define new variables (j)~2 = ^^^^ and tp^^ = where i = 1,3. As the outer junctions are 

identical, it is enough to study any one subsystem. So considering the case with i = 1, we 
write, 

'Pu = i'u (11) 

1 3 

i^i2 = -13^12 - COS sin + 2 [^dc + io cos{VLt)\ - as{-i)i2)- 

From Eq. [11] we conclude that in the presence of an external applied field it is not possible 
to synchronize all the three junctions due to the asymmetry induced by the applied fields. 
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However in the absence of an external field, an identical solution can exist for all the three 
junctions. Extending the symmetry analysis to a system of N JJs coupled in parallell with 
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^2 

FIG. 3: (a) and (b) shows that the junctions are phase correlated. (3 = 0.3, io = 1-2, u; = 0.6, i^c = 
0.3, = 0.37and 9 = O.Stt. 

nearest neighbour coupling, the second and the (A^ — 1)*'^ junction may have an identical 
solution for certain parameter values. Similarly, the third and the (A^ — 2)"°' junctions 
may have identical solutions and so on. Thus in the case of an array, from symmetry 
considerations we may deduce that N/2 solutions may exist if there are even number of 
junctions in the array and ^^^^ solutions will be present for odd number of junctions. The 
time series plot for an array of 7 and 8 junction is plotted in FiglH It can be observed 
from FiglD^a) that in an array of seven JJs the four solutions exists for the parameter range 
considered. The fourth junction has an independent solution. In FigjU^b) we have plotted 
the time series for 8 JJs. 

The presence of a phase difference between the applied fields changes the scenario com- 
pletely. On the apphcation of a small phase difference between the applied fields, the outer 
junctions desynchronize and all the three junctions are thus uncorrelated. But for sufficiently 
large values of phase differences, all the three junctions are found to be in phase synchro- 
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FIG. 4: (a)shows the time series plot for an array of 7 JJs and (b) for 8 junctions. (3 = 0.3, Iq = 
1.2, uj = 0.6, idc = 0.3, as = 0.37and 9 = O.Svr. 

nization. Considering the difference variables '?/'i,2; "^1,3; '^'^3,2 as defined earlier, we explain 
the phenomena as follows. Due to the asymmetry that arises between the outer junctions in 
the presence of the phase difference we need the extra variable '?/^3,2 to analyze this situation. 
The equations for the three difference variables may be written by substituting eqE] as 

^12 = -P^i2 - COS </)^2 sin 0r2 + 2 I'^d^ + ^0 cos(^]t)] - ^si^ + ^^2) (12a) 

^13 = -(^^13 - COS sin + tQ sm{nt + ^) ~ "s(^r3 - ^32) (12b) 

1 'ip~ 
i^32 = -Pi^32 - COS (pt2 sin 032 + ^ t^'^'^ + ^0 cos(fit + 9)] - as{-^ + ^12), (12c) 

where = io sin |. Thus each subsystems experiences a different driving field with the same 



frequency but different phases. Due to the phase relationship between the driving fields, a 
definite phase relationship is found to exist between all three junctions. 

The level of mismatch of chaotic synchronization can be given quantitatively by taking 
the similarity function S{t) as a time averaged difference between the variables tpi taken 



with time shift r 



and 



mm mm 



1/2 



(13) 



(14) 



mim] msm'^" ' 

and searching for its minimum cr = minrS{T). If 4Ji{t) = ip3{t), then S{t) has a minimum 
value cr = for r = 0. If both ipi{t) and ip3{t) are independent then S{t) ^ 1 for all the 
time. Line 1 in Fig. [5] shows complete synchronization between the end junctions and line 
2 shows that the outer and middle junctions are desynchronized when no phase difference 
is present. A minimum of S{t) indicates the the existence of a time shift between the two 
variables related to the phase shift. The amplitudes are uncorrelated in this regime, but 
phase correlation is present as indicated by lines 3 and 4 in the presence of a phase difference 
between the applied fields. On the application of a phase difference of 7r/2 the dynamics 




0.5 



FIG. 5: (1) is S{t) for variables V'l&V's i (2) for ^i&^/;2- Line (1) shows complete synchronization 
while the other is not synchronized. Both these lines are plotted with no phase difference applied 
(3) and (4) gives the similarity function for the variables Vi&V's V'i<^V'2 in the presence of 
phase difference 6 = O.Svr. The second figure shows line 4 where the dip can be observed clearly. 

changes to periodic one as can be observed from Fig. [6] 
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FIG. 6: The variables corresponding to the three JJs are ploted against time which indicates the 
periodic behavior. 

IV. RESULT AND DISCUSSION 

We consider a parallel array of JJs with linking resistor Rs and the conditions for synchro- 
nization is discussed. The outer junctions being symmetric, can possess identical solution 
and hence may synchronize depending on the parameter values. Linear stability analysis is 
done to find the stability of the synchronous solution of the outer junctions. The sum of 
conditional Lyapunov exponents calculated for the outer sub-system is found to be negative 
indicating stable synchronous state. Prom symmetry considerations we show that all three 
junctions could be synchronized only in the absence of an external field. Similarly in an 
array of N Josephson junctions, N/2 identical solutions may exist if the number of junctions 
is even and solutions may exist if the number of junction is odd. In the presence of 
a small phase difference, the system desynchronizes due to the asymmetry induced by the 
phase difference. As the phase difference is increased, in the case of three junctions all the 
three junctions act as if they are driven by different driving fields having the same frequency, 
but different phases. A phase synchronization is observed between all the three junctions 
and the motion becomes periodic. Thus, suppression of chaos can be obtained in Josephson 
junction systems in the presence of a phase difference between the applied fields and this 
property may find applications in the working of devices constructed using JJs. 
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